The QCD phase diagram from Schwinger-Dyson Equations 
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We study the phase diagram of quantum chromodynamics (QCD). For this purpose we employ 
the Schwinger-Dyson equations (SDEs) technique and construct a truncation of the infinite tower 
of equations by demanding a matching with the lattice results for the quark-anti-quark condensate 
at finite temperature (T), for zero quark chemical potential (/x), that is, the region where lattice 
calculations are expected to provide reliable results. We compute the evolution of the phase diagram 
away from T = for increasing values of the chemical potential by following the evolution of the heat 
capacity as a function of T and /i. The behavior of this thermodynamic variable clearly demonstrates 
the existence of a cross-over for ^l less than a critical value. However, the heat capacity develops a 
singularity near /i ~ 0.22 GeV marking the onslaught of a first order phase transition characterized 
by the existence of a critical point. The critical line continues until /x ~ 0.53 GeV where Tc = and 
thus chiral symmetry is finally restored. 

PACS numbers: 25.75.Nq, ll.30.Rd, ll.15.Tk, 11.55.Hx 



I. INTRODUCTION 

The strong interaction sector of the standard model 
of particle physics involves a phase transition which is 
relevant to the evolution of the early universe. At low 
temperatures, the observable degrees of freedom of quan- 
tum chromodynamics (QCD) are color-singlet hadrons 
whereas at high energies or temperatures the interac- 
tion gets increasingly screened and hence weak, causing 
hadrons to break up into a new phase where the dom- 
inant degrees of freedom are the defining ingredients of 
perturbative QCD, namely, quarks and gluons. This is 
dubbed as confinement- deconfinement phase transition. 
The same physical picture prevails when chemical poten- 
tial is increased. Experiments as well as lattice QCD 
computations for physical quark masses suggest that the 
temperature driven transition at zero chemical potential 
jjL is not a thermodynamic singularity. Rather, it is a 
rapid but smooth crossover from the regime describable 
as a gas of hadrons, to the one characterized by quarks 
and gluons [l|. The ^ driven transition at zero tem- 
perature T is qualitatively different. This regime can- 
not be accessed through lattice studies because the com- 
putations are hampered by the notorious sign problem. 
Nevertheless, a number of different model approaches in- 
dicate that the transition in this region is strongly first 
order Q]. 

Chiral symmetry restoration is another phase tran- 
sition which is expected to occur when temperature 
and (or) chemical potential are sufficiently large. This 
happens when the dynamically generated component of 
quark masses vanishes. The pattern of chiral symme- 
try breaking at zero T and /i already provides evidence 
of the fact that large effective quark masses within light 
hadrons owe themselves primarily to the strength of the 
QCD interaction and dictate both their static and dy- 



namic properties, see for example a recent review 0]. 
Its experimental implications for the elastic and transi- 
tion form factors of mesons and baryons form an inte- 
gral part of the planned program at the 12 GeV upgrade 
of the Thomas Jefferson National Accelerator Facility in 
Virginia [J] . As the strength of the QCD interaction di- 
minishes with increasing T and /i, only the bare quark 
masses survive when these parameters exceed a critical 
set of values. This is referred to as the chiral symmetry 
breaking- chiral symmetry restoration phase transition. 

In connection with these phase transitions, there are 
key questions which, to a considerable extent, provide 
the motivation for the heavy ion collision program at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL), the Large Hadron Collider 
(LHC) at the European Centre for Nuclear Research 
(CERN) and the Compressed Baryonic Matter (CBM) 
experiment of the future Facility for Antiproton and Ion 
Research (FAIR) at Darmstadt. These can be summa- 
rized as follows : 

1. Since the first order line originating at T = in 
the QCD phase diagram cannot end at the /.t = 
axis which corresponds to the starting point of the 
cross-over line, it must terminate somewhere in the 
midst of the phase diagram. This point is gener- 
ally referred to as the critical end point. On one 
hand, it rests with the experiments to find observa- 
tional evidence for it through extracting chemical 
and thermal freeze out temperatures, and on the 
other, the value of a sound theoretical prediction 
for its existence and location on the QCD phase 
diagram can hardly be exaggerated. 

2. Are the two phase transitions, corresponding to de- 
confinement and chiral symmetry restoration, co- 
incidental? As they both owe themselves to the 



diminishing of the QCD interaction strength, one 
may expect them to chart out the same curve in 
the QCD phase diagram. However, its prediction 
and measurement continue to be essential as well 
as a challenging endeavor. 

In this article, we concentrate on the exploration of 
the first point. Monte Carlo methods of lattice com- 
putations are severely handicapped off the n = axis, 
because the fermion determinant becomes complex and 
thus standard Monte Carlo methods fail, as the integrand 
is no longer real and positive definite. However, these 
techniques can still be adapted to extract some informa- 
tion on the QCD phase diagram for /i ^ Q- On the 
other hand, the fundamental field theoretical equations of 
QCD, namely, Schwinger-Dyson Equations (SDEs), have 
no such restriction. As their functional derivation makes 
no appeal to the smallness of the interaction strength, 
they provide an ideal, natural and unified framework to 
explore all the nooks and corners of the QCD phase di- 
agram with ease. Moreover, the bare quark masses can 
be set as small as required and even studying the chirally 
symmetric Lagrangian is far from being herculean [6| . 

We study QCD phase diagram using the tools of SDEs 
at finite temperature and chemical potential. We work 
with two-fiavor QCD with physical up and down current 
quark masses assuming isospin symmetry. The simplest 
two-point quark propagator is a basic object to analyze 
dynamical chiral symmetry breaking and confinement. 
At finite T and /i, we start from the general form of the 
quark propagator 

S^^{p,Wn) = iA{p^,wl)j-p + i'j4wlC{p'^,wl) 

+ B{p\wl), (1) 

where Wn = Wn + ip- and u)„ ~ 2{n + 1)ttT are the Mat- 
subara frequencies. A{p^,wl^), B{p'^,w'^) and C{p^,w-^) 
are the scalar functions to be self consistently determined 
through solving the corresponding SDE 

S^'^{p,Wn) = ij ■p + ij4.Wn+^ip,Wn) , (2) 

where I](p, z«„) is the self-energy expressed in terms of 
the dressed gluon propagator D^^ip — q, ^ni), with Q.ni = 
Wn — wi, and the full quark-gluon vertex T^{q, wijp, Wn) 
as follows 



^ — — oo 

^ -7rlt^S{q,wi)—r^{q,wi,p,Wn). (3) 



Pj^^ being the transversal and longitudinal projectors 
defined as 



pL _ c _ kfjki, 



^U 1 
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(5) 



Following the lead of Qin et. al. Q, we define 
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D,iP,nl,) ^ ViT)^iP + ^l^)e-(i'+^l^+<)/^' (6) 

where cr = 0.5 GeV. There are a couple of noticeable 
differences in our adaptation of the above ansatz as will 
be clear along the way. Note that we take D^ 7^ Dt 
which is generally true when T, pi ^ 0. The difference 
is governed by the gluon Debyc mass nig which enters 
only the longitudinal part of the gluon propagator. We 
take its value to be the lowest order hard thermal loop 
resuh ml = 4g^T^/3 Q with 5 = 1. The effect of the 
Debye mass is negligible all the way through the ther- 
mal evolution, even when we use larger values of g. The 
role of the Debye mass is to help for the convergence of 
the solutions. We also allow I? to be a function of tem- 
perature. With appropriate traces and simplifications, 
we obtain three coupled integral equations for A{p^, ^n)j 



B{p'^,w'^) and C{p'^,w'^) which are solved numerically. 
Convergence tests for the results allow us to sum over 14 
Matsubara frequencies without sacrificing the accuracy 
of the results by more than 5%. Once we have calcu- 
lated the quark propagator, we can evaluate the quark 
anti-quark condensate as : 



W) 



N.TY^j 



d' 



q 



(27r) 



Tr[S{q,Wn)] 



(7) 



However, notice that we have not yet said anything ex- 
plicitly about V{T) (coming from the gluon propagator) 
or the quark-gluon vertex, or in fact their product which 
gets projected onto the quark self-energy. Where do we 
extract this information from? At this point, we intro- 
duce the most important input ingredient for our calcu- 
lations. As we expect lattice computations to provide 
reliable results along the /x = 0-axis of the QCD phase 
diagram, we make the following ansatz : 



V{T) T^{q,wi,p,w„) = D{T) 7^ 



(8) 



and deduce D{T) so as to reproduce the lattice results 
for the temperature dependent quark anti-quark conden- 
sate quoted in Q. Note that for the chiral limit, we 
reproduce the results presented in [7| with appropriate 
modifications. In order to introduce a finite bare mass 
for the quarks, we use the fact that the mass function 
can be written as M{p'^,w'^) = B{p'^ ,w'^)/A{p'^,w'D and 
g^D^,{k, ilni) = p1^^Dt{P ^^ni) + P^,DL{P,nni) ,(4) then make the replacement B{p^,wl) -^ B{p^,wl) + 



Here g^ is the QCD interaction strength and the gluon 
propagator has the following general form in the Landau 
gauge 




0.05 0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 0.275 

T[GeV] 




0.135 0.14 



0.145 0.15 0.155 

T[GeV] 



\ , I ^ 

0.16 0.165 0.17 



FIG. 1: Function D(r) = a-bT^ -ctanh{d-eT^) of Eq. dH) 
at /i = reproduces the lattice data Q] for the quark anti- 
quark condensate {^tp) very well, especially in the region of 
the phase transition. 



FIG. 2: Heat capacity -9t(i/''0) calculated from the SDE for 
D{T) given by Eq. ([9| as a function of temperature. Its max- 
imum at Tc ~ 0.153 GeV agrees well with [y|). 



A{p'^,w'f^)'m. The self consistent solution of the cou- 
pled equations subsequently yields the quark propagator. 
With a fairly simple form of the function D{T), i.e., 



D{T) 



bT^ - ctanh{d - eT^) 



(9) 



at /i = (a = 2.17,6 = 343.64, c = 1.76, d = 0.78 and 
e = 273.8 with appropriate mass dimensions in GeV), 
lattice data is sufficiently well reproduced as shown in 
Fig. [T] One can draw the corresponding heat capacity 
curve given by — 9t('0'0) which yields a critical temper- 
ature of Tc « 0.153 GeV, see Fig. [5] (compare it with 
Tc = 0.154 GeV reported in @). 

We now extend our model and employ the same D{T) 
as we evolve our results for /x 7^ 0. The gluon Debye mass 
becomes ^-dependent, i.e., m^ = {4T^ /3 + ^i'^ /TT^)g^ , also 
with g = 1. For increasing /.t, we naturally observe chi- 
ral symmetry restoration at lower temperatures. This is 
depicted in Fig. |3] which confirms that the effect of in- 
corporating the chemical potential /U, tends to restore the 
chiral symmetry for lower temperatures. An important 
point to note is that as /x gets larger, a discontinuity 
starts to set in for the region where the critical tempera- 
ture is located. It becomes increasingly pronounced with 
a relatively small variation of /j,. Moreover, the slope of 
the curve increases drastically, so much that we even- 
tually observe a vertical drop, as is evident in Fig. |31 
This is a likely indicator of a physical effect. In order to 
quantify this change in the behavior of the heat capac- 
ity, we draw it as a function of temperature for different 
values of /i to locate the critical point. The result is 
shown in Fig. |4l The plots become narrower for increas- 
ing fi which helps locate the critical temperature with 
progressive precision. On the other hand, the height of 
the curve soon shoots up to infinity as evident from Fig. 2] 
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FIG. 3: The condensate for different values of /^ as a function 
of temperature. For large n, the curve develops a disconti- 
nuity which becomes more and more marked for increasing 
values of fi. 



before the chemical potential reaches a value of about 
0.22 GeV. We identify this thermodynamic singularity 
with the onslaught of a first order phase transition. Be- 
fore that point is reached, the phase transition is a cross- 
over. Notice from Fig. [3] that chiral symmetry remains 
broken beyond fj, « 0.22GeV, where the first order phase 
transition is triggered. Our numerical accuracy allows us 
to explore the chemical potential up to about 0.3 GeV. 
However, the trajectory of the points along the critical 
line is smooth enough to make a quadratic fit to estimate 
the value fi w 0.53 GeV for zero temperature where chi- 
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FIG. 4: We plot the heat capacity curve for the condensate as 
a function of temperature for different values of fj,. The peak 
gives the critical point {fic,Tc). Note that the height of this 
thermodynamic variable shoots up to infinity for a sufficiently 
large /i, indicating a change in the order of phase transition. 



ral symmetry is finally restored, as shown in Fig. [5] All 
these characteristics are sketched in Fig. El 

A few words of caution are here in order. The func- 
tion D{T) in Eq. (|S]) and the full quark-gluon vertex 
T^{q, wi^p, Wn) surely depend on /i and a different choice 
of these quantities may change the position of the criti- 



FIG. 6: This plot shows different regions of the QCD phase 
diagram as a function of the chemical potential n. From ^ = 
to the vertical dashed line at /i « 0.22 GeV, the transition is 
a cross-over. The height of the heat capacity, the thermody- 
namic variable drawn on the vertical line, becomes singular 
as it approaches the dashed line, suggesting a change in the 
nature of phase transition from a simple cross-over to the first 
order phase transition. The chiral symmetry continues to be 
broken above /i ~ 0.22 GeV, as evidenced from Fig. (3] How- 
ever, it becomes numerically demanding to find the transition 
points. Chiral symmetry appears to be restored at ^ ~ 0.53 
GeV, represented by a vertical dot-dashed line. 
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FIG. 5: We plot the phase diagram of QCD, indicating also 
the critical end point which corresponds to fj. ~ 0.22GeV and 
T = 0.097GeV . Due to numerical difficulties above around 
fi ~ O.SGeV, we adopted a quadratic numerical fit /i + f2T + 
fsT^, to extrapolate our results to the T = axis of the 
phase diagram which yields /i ~ 0.53 GeV above which chiral 
symmetry is restored. 



cal end point. Moreover, the fermion-boson vertex must 
be constructed from theoretical [lOl, as well as phe- 
nomenological [ll| constraints. At finite T and (or) fi, 
the complexities involved are highly non-trivial but work 
has already begun in this direction |12| . We have per- 
formed preliminary studies to explore the nature of the 
confincmcnt-deconfinement transition and have found 
that within this model the transition line follows the foot- 
steps of the chiral symmetry breaking-restoration transi- 
tion. This indicates that other transitions to quarkionic 
matter or a color locked phase, as the chemical potential 
increases, can only be described by incorporating addi- 
tional ingredients to this model. While the above men- 
tioned details will certainly be trimmed over the next few 
years, this letter shows that SDKs arc an efficacious tool 
to explore the QCD phase diagram. We have used this 
tool to locate the critical end point as well as to pro- 
vide clear signals of the cross over (near ^ = 0), and of 
a first order phase transition as we move away to suffi- 
ciently high values of /i, starting from lattice results for 
the finite temperature quark-antic^uak condensate. 
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